Stopping in central Pb + Pb collisions at SPS energies and beyond 



Yacine Mehtar-Tani 1 and Georg Wolschin 2 

1 Departamento de Fisica de Particulas and IGFAE, Universidade de Santiago de Compostela, Spain 

2 Institut fur Theoretische Physik der Universitdt Heidelberg, Philosophenweg 16, D-69120 Heidelberg, Germany 



o 

C 



PACS 25 . 75 . -q - Relativistic heavy-ion collisions 

PACS 24.85.+p - Quarks, gluons, and QCD in nuclear reactions 

PACS 25.75.Nq- Quark deconfmement, quark-gluon plasma production, and phase transitions 

Abstract. - We investigate stopping and baryon transport in central relativistic Pb + Pb and Au 
+ Au collisions. At energies reached at the CERN Super Proton Synchrotron (tJsnn = 6.3-17.3 
GeV) and at RHIC (62.4 GeV), we determine the fragmentation-peak positions from the data. 
The resulting linear growth of the peak positions with beam rapidity is in agreement with our 
results from a QCD-based approach that accounts for gluon saturation. No discontinuities in the 
net-proton fragmentation peak positions occur in the expected transition region from partons to 
hadrons at 6-10 GeV. 
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Introduction. — The search for signatures of the 
, <~| cross-over of the nuclear medium at high temperature and 
l— — low chemical potential from a confined to a deconfined 
(-s^ phase in relativistic heavy-ion collisions as predicted by 
lattice QCD [Ij, and for a first-oder phase transition at 
■^j" lower temperatures - with a critical point at intermediate 
temperatures - has a varied history. Some of the decisive 
^^observables in this context such as strangeness enhance- 
• ment [2] have been discussed since more than 20 years. 

f***) Energy scans of such observables are indeed suitable 
t— (to try to detect the threshold energy for the QCD phase 
transition. As an early example, a maximum in the 
j>< K + > / < 7T+ > ratio was found by the NA49 col- 
laboration in central Pb + Pb collisions at Jsnn — 8 



GeV which is "... consistent with the hypothesis that a 



X 

C3 transient state of deconfined matter is created above these 
energies" 3). There exist, however, other interpretations 
of these data. In particular, the maximum is consistent 
with the decreasing net-baryon value combined with in- 
creasing KK pair production at rising ^s^n [Ij. 

The ongoing STAR beam energy scan program BES at 
RHIC [5] tests the low-energy range y / s]v]v = 5 — 50 GeV 
in Au + Au collisions with the aim to look for specific 
signatures of a critical point in fluctuations and collec- 
tive effects using data for identified particles from a single 
detector. It is also expected to see the disappearance of 
partonic effects at lower energies. 

For </>-Meson production, data from NA49 and STAR 
exist already for central Cu + Cu, Au + Au and Pb + Pb 



collisions in the energy range 6.3 — 200 GeV. It was argued 
in [6] that these data "indicate a threshold energy for the 
confinement-deconfinement phase transition" at ^/sjvw = 
15.7 ±8.1 GeV. 

It has also been proposed to look for signatures of the 
deconfmement threshold in the net-baryon (B — B) or net- 
proton rapidity distributions. As opposed to distributions 
of produced charged hadrons which are close to statistical 
equilibrium, these are clearly nonequilibrium distributions 
with two fragmentation peaks in rapidity space that move 
away from midrapidity with increasing cm. energy [7j[8]. 
For a thorough test of theoretical models one needs the 
full rapidity distribution, not just the midrapidity region. 

In (9j[l0) it was shown on the basis of a Relativistic 
Diffusion Model 18] that the net-proton rapidity distribu- 
tions in central Au + Au at RHIC energies of 200 GeV 
can not be understood as being due to hadronic diffusion 
processes, and that the midrapidity region is populated 
through partonic processes that indicate deconfmement. 
In contrast, at AGS and low SPS energies the measured 
net-proton rapidity distribution can be modeled precisely 
within the hadronic diffusion description including a pre- 
diction of the correct energy dependence. 

The energy dependence of the net-proton rapidity distri- 
butions at AGS and SPS energies has also been discussed 
by [ll] with respect to signs for a transition to the quark- 
gluon phase: A possible non-monotonous behaviour of the 
curvature of this distribution at midrapidity as function of 
the incident energy was investigated. There is an indica- 
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tion for such a behaviour in the region ^snn — 4-10 GeV 
from three-fluid dynamic calculations with an equation of 
state that involves a first-order transition 11 , but it is 



not clear or even doubtful whether the alleged curvature 
anomaly is supported by the data with sufficient statistical 
significance. 

In this Letter we investigate the net-proton distributions 
for central Pb + Pb collisions in the CERN Super Proton 
Synchrotron energy region ^/satjv = 6.3-17.3 GeV, and 
central Au + Au at RHIC energies of 62.4 GeV. 

We have made a detailed analysis of stopping in Pb 
+ Pb collsions at 17.3 GeV with NA 49 data from 
12 13 , and of Au + Au collisions at RHIC energies with 
BRAHMS data [IJJ already in 15fl6 . Our partonic model 
accounts for interactions of the valence quarks with the 
gluon condensate in the respective other nucleus. It is not 
meant to fully account for stopping at considerably lower 
energies. The aim is here to look for indications of the 
hadron-parton transition in net-proton rapidity distribu- 
tions at SPS energies through deviations between model 
results and data at low energies, or discontinuities as func- 
tion of energy. 

With this aim in mind, we use the fragmentation-peak 
position as function of beam rapidity, or center of mass 
energy, as an observable. Since it can be determined quite 
accurately from the data at low energies, and is expected 
to depend linearly on the beam rapidity from the partonic 
model, it is likely to be a more precise indicator than the 
mean rapidity loss that had been discussed in 17 . The 



mean rapidity loss deviates from a linear dependence on 
beam rapidity at RHIC energies and beyond, in accor- 



dance with our partonic model 15 16 



We reconsider essential ingredients of the model in 
the next section, where the linear dependence of the 
fragmentation-peak positions on beam rapidity that the 
model predicts is emphasized. Then we determine the 
fragmentation-peak positions in analyses of the SPS and 
RHIC data, and find good agreement with the model pre- 
dictions, but no evidence for discontinuities. We draw the 
conclusions in the last section. 

Stopping in QCD. — The differential cross section 
for valence quark production in a high-energy nucleus- 
nucleus collision is calculated from 
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dN 



1 1 



<P PT dy (2tt) 2 p| 



xiq v (xi,Qf) <p(x2,pr), (1) 



where px is the transverse momentum of the produced 
quark, and y its rapidity. The longitudinal momen- 
tum fractions carried, respectively, by the valence quark 
in the projectile and the soft gluon in the target are 
xi — pr/*/sexp(y) and X2 — pr / ^ exp(-y) . The fac- 
torization scale is set equal to the transverse momentum, 
Qf = pt- The valence-quark distribution is q v (xi,Qf). 
We have discussed the gluon distribution (p(x,pr) and de- 
tails of the overall model in 15 16 . 
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Fig. 1: Double-gaussian fits of the net-proton rapidity distri- 
butions in Pb + Pb (from top to bottom) at ^/smn = 6.3, 
12.3 and 17.3 GeV, and in Au + Au at 62.4 GeV allow to de- 
termine the positions of the fragmentation peaks in rapidity 
space. Preliminary Pb data points are from the NA49 collab- 
oration |13||18] using a conversion factor 82/208 from baryons 
to protons. At 17.3 GeV the data shown here are from the 
previous NA49-analysis [12] . Au data in the bottom panel are 
from BRAHMS [531. Error bars are smaller than the point 
sizes unless shown. See table [T] for results at all investigated 
energies, including the new Pb data at 17.3 GeV. 
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Fig. 2: Peak positions of the net-proton rapidity distributions 
in Pb + Pb (black squares) and Au + Au (black diamond) 
as function of the beam rapidity, determined from double- 
gaussian fits of the NA49 
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see ta- 



13||18] and RHIC data 
ble [I] The open square is based on older NA49 data at 17.3 
GeV [12]. The slope agrees well with the analytical expression 
eq.|3| for A = 0.2, dashed line. The cross refers to the cal- 
culated peak position of Au + Au at 200 GeV, the inclined 
crosses to Pb + Pb at LHC energies of 2.76 and 5.52 TeV. 



One important prediction of the color glass condensate 



theory ,22-25 is geometric scaling: the gluon distribution 
depends on x and pt only through the scaling variable 
Pt/Q^( x )> wn ere Q 2 s (x) = A 1 / 3 Ql x~ x , A is the mass 
number, Qq sets the dimension, and A is the saturation- 
scale exponent. This has been confirmed experimentally 
for A = 1 in deep inelastic e+p scattering at HERA 26 . 
There the fit value A = 0.2 — 0.3 agrees with theoretical 
estimates based on the Balitsky-Kovchegov equation (27l 



28 that include a running coupling 29-31 . 

To show that the net-baryon distribution in relativistic 
heavy-ion collisions also reflects the geometric scaling of 
the gluon distribution, the following change of variables 
has been proposed in 
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-2y 



Pt 



x = Xi, x-i = x e 
x's e~ 2y . Thus after integrating over px, the rapidity dis- 
tribution can be rewritten as 



dN 
dy 



(r) 



C 

2n 



1 dx 



xq v (x) if(x + e T ) 



(2) 



where r = ln(s/<3o) — In A 1 / 3 — 2(1+A) y is the correspond- 
ing scaling variable and C is the overall normalization that 
contains the fragmentation function |16| , and is fitted here 
to the experimental yields. Hence, the net-baryon rapidity 
density in the peak region is a function of a single scal- 
ing variable r, which relates the energy dependence to the 
rapidity and mass number dependence. 

For the fragmentation-peak position, we have derived in 
15 the analytical result 



Vpeak 



1 

r+A 



lib, 



h\A 1/G ) + const. 



(3) 



Table 1: Energies and beam rapidities for Pb + Pb (Au 
+ Au a ), and corresponding positions and widths of the net- 
proton fragmentation peaks from double-gaussian fits. The fit 
results Vpea k an< ^ ^peak (FWHM) with estimated error bars re- 
fer to the NA49 data [l3|[T8] and to the BRAHMS data 
Results marked as 



are from a fit to previous NA49 data 



(GeV) 


\ljbeam \ 


\ypeak 1 


p/it 
peak 


6.3 


1.92 


0.67±0.10 


1.78±0.2 


7.6 


2.10 


0.80±0.10 


2.08±0.2 


8.8 


2.25 


0.91±0.25 


2.19±0.5 


12.3 


2.58 


1.14±0.15 


2.25±0.3 


17.3 


2.92 


2.00±0.40 


3.11±0.6 


17.3 6 


2.92 


1.26±0.25 


2.17±0.4 


62.4 a 


4.20 


2.63±0.10 


3.33±0.2 



It is expected to be applicable at RHIC energies and 
beyond, but here we test the validity of the predicted linear 
growth with y\, eam at SPS energies. At very low energies 
deviations will occur due to kinematical constraints in the 
limit ybeam ~^ where y pe ak -> 0. 

Results. — Using the model as outlined in the pre- 
vious section, we have investigated baryon transport in 
central Pb + Pb collisions at SPS energies of -Jsnn — 
6.3-17.3 GeV, and in central Au + Au at RHIC ener- 
gies of 62.4 and 200 GeV. In the high-energy region, our 
parton-based approach gives good results for net-proton 



rapidity distributions 15 16 with respect to the available 
BRAHMS data [7)[l4]. 

With the parameters fixed at 200 GeV (A = 0.2, and 
Qq — 0.04 GeV 2 ), we have also used it for predictions at 

These theoret- 
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LHC energies of 2.76 and 5.52 TeV 
ical results at LHC energies will, however, not be tested 
experimentally in the next years at or near the fragmen- 
tation peaks since particle identification at these energies 
is restricted so far to \y\ < 2. 

In the region of SPS and low RHIC energies, we have 
determined the net-proton fragmentation-peak positions 
Vpeak from the preliminary data of the NA49 collabora- 
tion [13[[18] and of BRAHMS JH] in double-gaussian fits, 
as shown in fig. [l] and table There is an uncertainty 
regarding the peak position at 17.3 GeV since the experi- 
ments from two different runs at this energy have produced 



different results 12 18 , with the more recent preliminary 
data not extending across the fragmentation-peak posi- 
tion. 

Different from the Au + Au data at 62.4 GeV, the 
experimental BRAHMS results at 200 GeV [7] can not 
be reproduced in a double-gaussian fit to determine the 



p-3 



Y. Mehtar-Tani and G. Wolschin 




Fig. 3: Widths (FWHM) of the net-proton fragmentation-peak 
rapidity distributions in Pb + Pb (squares) and Au + Au (di- 
amond) as functions of the beam rapidity based on double- 
gaussian fits (fig. flj) of the available SPS data (black squares) 
from NA49 [12|13|18] , and of RHIC data from BRAHMS fffli] . 
The widths increase at SPS energies. At high RHIC energies 
and at LHC, a constant value is expected from our model. 



fragmentation-peak positions. Instead, an additional 
midrapidity source for net protons appears that had been 
interpreted in [9j as being due to partonic contributions 
to the net-proton yield. It is related to valence quark 
- gluon collsions since gluon-gluon collisions do not con- 
tribute to net-baryon distributions. Indeed the data are 



in good agreement with our partonic calculation 15 16 
which had been calibrated at this energy. 

The peak positions determined from the fits as func- 
tions of beam rapidity are shown in fig. [2] Here the beam 
rapidity is obtained from the center-of-mass energy as 
ybeam = hi^ s n n / rn p ) with the proton mass Trip . The 
peak positions are found to increase linearly with beam 
rapidity at SPS and RHIC energies. The results of the 
double-gaussian fits are in surprisingly good agreement 
with the prediction of our microscopic model in eq. Q 
for a saturation-scale exponent A = 0.2, an empirical con- 
stant = -0.2, and A=208 (208 1 / 6 ~ 197 1 / 6 ). 

At SPS and RHIC energies, a saturation-scale exponent 
A = 0.2 clearly gives the best result for the slope in com- 
parison with the NA49 and BRAHMS data, in agreement 
with results of [3l]|33] . At high RHIC energies of 200 GeV 
and LHC energies of 2.76 and 5.52 TeV which are also 
indicated in fig. [2j a larger value of A may turn out to be 
more adequate, depending on the growth of the saturation 
scale with energy. 

A saturation-scale exponent A = 0.2 corresponds to a 
gluon saturation momentum of Q s ~ 0.77 GeV at a;=0.01. 
When comparing this value to investigations of charged- 



hadron production such as in 31 , which involve the gluon 



corresponding to a factor 9/4 [32]. This leads to Q s ~ 1.2 
GeV at x=0.01. 

There is no discontinuity in the peak positions as func- 
tions of center of mass energy, nor any disagreement be- 
tween calculation and data analysis that might indicate 
the sudden onset of purely hadronic processes at lower 
energies. Hence, the transition from partonic to hadronic 
processes in the observable net-proton fragmentation-peak 
positions is gradual, without immediate signatures of a 
cross-over as predicted from lattice gauge theories for equi- 
librium systems of partons. 

The widths of the net-proton distribution functions as 
determined from the SPS and RHIC data in the double- 
gaussian fits are shown in table [T] and fig. [3] Plotted are 
the full widths T at half maximum (FWHM) of the indi- 
vidual fragmentation peaks, which are related to the stan- 
dard deviations a through T — \/81n2 • a. Similar to the 
mean values, they increase with beam rapidity from SPS 
to RHIC energies. The widths will tend to the width of 
the Fermi distribution in rapidity space as ybeam — > 0. 

An increase of the width with center-of-mass energy is 
expected in the Relativistic Diffusion Model [34] that has 
been used to interpret SPS and 62 GeV RHIC data. In 
the high-energy limit, our microscopic model predicts the 
width to stay constant, but again this result will be dif- 
ficult to test at LHC energies. The transition from an 
increasing width at SPS energies, to a constant width at 
higher RHIC energies and beyond, would be indicative of 
the transition from predominantly hadronic to partonic 
system behaviour. 

Conclusion. — An analysis of net-proton rapidity dis- 
tributions in central Pb + Pb collisions at SPS cm. ener- 
gies of 6-17 GeV, and Au + Au at low RHIC energies of 
62.4 GeV has been presented. The fragmentation-peak po- 
sitions have been determined in double-gaussian fits of the 
available data. On the theoretical side, the investigation 
uses our microscopic approach based on quantum chro- 
modynamics that analytically predicts the fragmentation- 
peak positions as function of center-of-mass energy, or 
beam rapidity 



15 16 



distribution in the adjoint representation of SU(3), one has 
to consider a rescaling of our net-baryon Q 2 S by the color 
factor N C /C F with C F = (iVg - l)/2N c and N c = 3, 



When tested at low energies, good agreement with the 
fragmentation-peak positions determined from the data in 
the measured net-proton distributions occurs at cm. ener- 
gies of 6-62 GeV. The linear growth of the fragmentation- 
peak position with beam rapidity that our microscopic 
model predicts is indeed found in the SPS and RHIC data. 

There is no discontinuity as function of energy, nor a 
severe disagreement with the data that might have indi- 
cated the sudden onset of hadronic processes at low ener- 
gies as consequence of a parton-hadron transition. Hence, 
the highly non-equilibrium net-baryon peak positions in 
rapidity space do not seem to be well-suited as signatures 
for the parton-hadron cross-over, or first-order phase tran- 
sition at higher values of the chemical potential. 

There are indications for a cross-over signature in net- 
proton distributions near midrapidity, as has been sug- 
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11 . However, there are probably better [30] 



gested in 

chances to establish the existence of the hadron-parton 
transition in obscrvables that are much closer to statisti- 
cal equilibrium, such as rapidity distributions and cross- 
sections of produced strange mesons. The fact that such 
observables approach statistical equilibrium much closer 
than net-baryon distributions makes it more plausible to 
compare with lattice results regarding the hadron-parton 
cross-over. 
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